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Abstract 

Spatiotemporal chaos (STC) exhibited by the Kuramoto-Sivashinsky (KS) 
equation is investigated analytically and numerically. An effective stochastic 
equation belonging to the KPZ universality class is constructed by incor- 
porating the chaotic dynamics of the small KS system in a coarse-graining 
procedure. The bare parameters of the effective theory are computed ap- 
proximately. Stability of the system is shown to be mediated by space-time 
defects that are accompanied by stochasticity. The method of analysis and 
the mechanism of stability may be relevant to a class of STC problems. 
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I. INTRODUCTION 



A major goal in the study of spatiotemporal chaos (STC), the chaotic dynamics of a 
spatially extended deterministic system, is to obtain a statistical description in the limit 
of large length and long time scales (the hydrodynamic limit) [1,2]. One avenue towards 
this end is the use of effective stochastic equations. Once effective equations are established, 
standard methods of statistical field theory can be employed to described the behavior in the 
hydrodynamic limit [3]. Some outstanding issues include understanding when the stochastic 
descriptions are applicable and how hydrodynamic equations may be constructed. 

Here, we construct an effective hydrodynamic theory for the Kuramoto-Sivashinsky (KS) 
equation in 1 + 1 dimensions [4-6], 

d,h = -vdlh - Kd^h + ^{d,h)\ (1.1) 

which, due to its simplicity, has often been regarded as a paradigm for STC. The KS equation 
(with ;/,/€> 0) is believed to describe a variety of systems driven away from equilibrium. 
Among them are fiame front propagation [4], the fiuctuations of a fiuid fiowing down an 
inclined wall [5], and phase turbulence in convection cells [6]. 

The qualitative phenomenology of STC in the KS equation is rather simple: The system is 
linearly unstable to long-wavelength perturbations, and locally forms "cells" of a preferred 
size, given by the most unstable wavelength. It is however (exponentially) unlikely for 
these cells to be in complete registry in a large system [7] because they initially develop 
independently in different regions according to the linearized equation. Generically, the 
cells are locally compressed or stretched from the preferred size. Too much compression 
eventually leads to cell annihilation in the nonlinear regime; similarly too much stretching 
leads to cell creation [2]. These annihilation and creation events are "space-time defects" 
formed to reduce the local cell misalignment. However, the defects generally do not improve 
cell alignment at the global scale but serve to generate more defects. Thus the chaos is 
characterized by local cellular patterns interspersed with dislocations [8]. Figure 1 shows 
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a spacetime profile of tlie cliaotic state. Defect-mediated or pliase turbulence is a term 
sometimes used to describe this situation [1]. 

Previous studies of the KS equation include aspects of the statistical properties [9-12], 
the spectrum of Lyapunov exponents [13], the attractor /bifurcation map for small sys- 
tems [14,15], and the stability analysis of the stationary solutions [7,16]. Ref. [1] provides 
a recent review. It was conjectured by Yakhot [17] that in the hydrodynamic limit, the 
probability distribution of the KS equation may be mimicked by that of the noise-driven 
Burgers' equation [18], 



The effective parameters z/gff and Dgff are assumed to be positive constants, and Agfr = A. It 
is useful to think of h as describing the height profile of a driven interface. In this context, 
Eq. (1.2) is also called the KPZ equation [19,20], whose solution's are well known in one 
dimension. 

Yakhot's conjecture was motivated partly from symmetry considerations. The KPZ 
equation is the simplest equation which preserves two important symmetries of the KS 
equation, the translational symmetry h — )• Const, and an infinitesimal Galilean symmetry 
h ^ h-\- 6x, X ^ X -\- Agfr^t, with Agfj = A. In addition, the apparent stability of the solution 
at large scales demands a positive surface tension, z/gfj , while the chaotic degrees of freedom 
motivates a Langevin noise. This conjecture was found to be consistent with results of large 
scale numerical simulations of the KS equation [11,12]. However, an explicit demonstration of 
the connection between KS and KPZ has been lacking. In particular, the precise microscopic 
mechanism of generating stability (z/gfj > 0) is not well understood. 

In this paper, we show that the stability of the KS system at large scales is mediated 
by the space-time defects. We have dubbed this phenomenon "defect-mediated stability". 





where rj is an uncorrelated Gaussian noise, with second moment 
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The chaotic nature of the defect dynamics then imphes that stabihty and stochasticity are 
unavoidably interconnected. To obtain these conclusions, we constructed an effective hydro- 
dynamic theory by explicitly coarse-graining the KS equation and utilizing the statistical 
properties of the chaotic dynamics in small KS systems (on the order of a few cell sizes). 
The method and the stability mechanism may be relevant to a class of STC problems. 

We begin in Sec. II with a review of some basic properties of the KS equation. In Sec. Ill, 
we re-examine the perturbation theory about the stationary cellular solutions, where space- 
time defects are excluded by definition. We illustrate how the dynamics of a large system 
may be obtained from various properties of small systems, and show the lack of diffusive 
stability (surface tension) in the large scale dynamics to be the direct consequence of a 
conservation law and symmetries pertaining to the cellular state. In Sec. IV, we extend the 
analysis to the chaotic regime and show how a stabilizing effective surface tension z/gff > is 
generated from the occurrence of space-time defects. We propose a coarse-graining procedure 
which allows us to construct the effective hydrodynamic equation from the chaotic dynamics 
of small KS systems subject to external forcings. Relevant statistical properties of the small 
system are obtained numerically in Sec. V. In Sec. VI, we analyze the results which lead 
to the KPZ equation as the effective stochastic equation in the hydrodynamic limit, and 
yield the approximate values of the bare KPZ parameters. A summary of the main ideas, 
the methods, the key assumptions, and the major results is given in Sec. VII, where we also 
discuss possible ramifications of our approach to other STC problems. A review of stochastic 
dynamics is given in the Appendix. 



We begin with a review of the linear stability analysis. For small fiuctuations in /i, the 
nonlinear term in Eq. (I.I) can be neglected. In this case, the linear equation exponentially 



amplifies all modes with wavenumber k < ki, where ki = Jv/ k. The fastest growing mode 



is kQ = a/z//2/€, corresponding to a preferred wavelength of Iq = 2tt jkQ. It has a growth time 



II. BASIC PROPERTIES 
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to = Ak/v'^. Thus, small fluctuations in the initial conditions become of order one on the 
order of the linear growth time to. Beyond this time the nonlinear term becomes important. 
It causes a DC "drift" in the solution since dth = ^[d^hy ^ 0. Here the overbar denotes 
spatial average over the entire system with periodic boundary conditions. Without loss of 
generality, we consider A > throughout the text. The time averaged drift rate, 

VL = (TO (2.1) 

is then a finite positive constant for a system of size Z, i.e. the solution drifts "upward" to 
larger values of h. The large scale structure, or roughness, of the interface profile originates 
from fiuctuations in the /oca/ drift rate about its average, v^. Note that the three parameters 
in Eq. (1.1) define the basic space, time, and height scales, /o,to, and ho = 2p / X respectively. 
The length of the system L/lo is the only dimensionless parameter. In the hydrodynamic 
limit L/lo oo, the KS equation can be considered parameter free. 

Consider first small systems with periodic boundary conditions where the drift rate is 
spatially uniform so that the dynamics are simple. For a system of size L < li ^ 27r/A;i, 
there are no unstable modes, so the trivial solution h = Const is a global attractor. For 
Z ^ /i, a nontrivial stationary, cellular solution H{x;L) appears, with L parameterizing 
the dependence of the solution on the system size. The cellular solution is defined by the 
ordinary differential equation, 

in a periodic box of size L. The drift rate Wo{L) is explicitly space and time independent for 
these solutions. It is found [14,16] that for li < L < I2 ^ 1.27/o, Fq. (2.2) admits a unique 
uni-cellular solution which is the global attractor. For L > multi-cellular solutions are 
possible. The KS equation undergoes a complicated bifurcation sequence with regions of 
stable multi-cellular solutions between regions of chaos as L is increased beyond I2 [14,15]. 
The basins of attractions of the stable solutions become vanishingly small for L ^ G/q. 

The stability of the stationary cellular solutions in the region li < L < I2 have been 
investigated in Ref. [16]. (Ref. [16] actually studied the stability of the velocity field u = d^h^ 
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but the results carry over directly to the potential field h considered here.) The precise 
shape of the cells depend on L/Iq. Figure 2 shows some examples of the cellular solution. 
For Z /q, one finds that smaller cells have larger amplitudes while larger cells have smaller 
amplitudes. Since the drift rate 



is a measure of the "non-fiatness" of the interface, the dependence of cell shape on L/Iq gives 
an increasing Wo[L) for decreasing L. On the other hand, the cellular solution must vanish 
at L = /i, so Wo[L) must decrease as L approaches li. The behavior of Wo[L) is shown in 
Fig. 3. It has a maximum at Z = /* 0.85/o. The change in drift rate upon a change of cell 
size is 



which changes sign at L = I*. As we will show, a also governs the response of the local drift 
rate to cell compression in a large system, and is crucial in determining the "elasticity" and 
stability of the large system. 

The existence of a continuum of stable, stationary uni-cellular solutions has profound 
effects on the properties of larger systems. For Z > 2/i, stationary cellular solutions com- 
prising of a periodic array of uni-cellular solutions H[x; I) can exist provided that the system 
size L is commensurate with the cell size /. This is not much of a constraint when there exist 
a broad continuum of Vs. Consequently, a degeneracy of periodic cellular solutions can exist 
for systems of lengths as small as a few /q's. For example, an L ^ 4/o system can fit four 
cells with size / /q or five cells of size / sa O.S/q. If we evolve the KS equation starting from 
random initial conditions, the system becomes frustrated attempting to choose between the 
four-cell or five-cell state. Such "frustration" effects can lead to chaotic dynamics [21], since 
it involves nonlinear coupling of at least three degrees of freedom: the amplitudes of the 
four-cell state, five-cell state, and their relative phase. The chaotic oscillations between the 
four-cell and five-cell state then lead to the "creation" and "annihilation" of cells, and are 
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perceived as "space-time defects" (Fig. (IV). Thus the dynamics (or statistical properties) 
of the defects can be obtained through an investigation of the small chaotic systems. This 
may be carried out using amplitude equations [21,22] or by direct numerical simulations of 
small systems. Here we shall assume that the relevant statistical properties of the small 
chaotic systems are already known, and derive the effective equation of motion describing 
the KS dynamics for L ^ /q. 

III. PERTURBATIONS ABOUT THE CELLULAR SOLUTIONS 

The effective hydrodynamics describing perturbations of the cellular solutions are derived 
first. They have been investigated by Frisch, She, and Thual [16] formally using a multiple- 
scale perturbation analysis. A similar study was done by Shraiman [7]. Here we shall 
recover these results from a simple coarse-graining procedure which highlights the underlying 
physics. Key concepts and methods developed here will be carried over directly into the 
chaotic regime in Sec. IV. Appreciating the physics of the cellular state will be crucial 
towards developing a theory for the chaotic state. 

We consider the effect of long-wavelength perturbations of the cellular solution H[x;£) 
with period li < £ < I2 so that the cellular solution is locally stable. The /i-field is divided 
into a "slow" part and a "fast" part /i>, where slow and fast correspond to long and 
short spatial distances with respect to £. More generally, we define 

[OU{x,t) = j£^'dyO{y,t), (3.1) 
and [0]y{x,t) = 0{x,t) - 0<{x,t), (3.2) 

where the operator (9 is a functional of h. We shall be interested in the long wavelength 
fiuctuations in /i<, for which 

dm<^[d:hU and [/.>]< ^ 0. (3.3) 

Spatially averaging Fq. (1.1) over £ and applying the approximations stated in Fq. (3.3), we 
obtain the equation of motion for the slow modes, 
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d,h^ = -vdlh^ - Kdth^ + ^{d^h^f + w^{x,t), (3.4) 



where 



w4x,t)^-J^ -fidyhyY (3.5) 

gives the coarse-grained local drift rate. 

To complete the coarse-grained dynamics, we need to obtain the equation of motion for 
w^. This can be derived formally by combining Eq. (3.5) with the equation of motion for 
the fast variable /i>, 

dthy = -vdlhy - nd^hy ^{d^hy fy Xd^h^d^hy. (3.6) 

Equation (3.6) is obtained by subtracting Eq. (3.4) from Eq. (1.1) since h = h^-\- hy. 

However, to emphasize some essential concepts that will be useful later we provide a 
physically motivated derivation. Consider fast modes comprised of slight local distortions 
of a cellular solution H{x] £), where £ parameterizes the shape and size of the solution. This 
is implemented through the use of the ansatz [23] 

hy{x,t) = H{x-l{x,t)), (3.7) 

where l[x^t) is the local cellular wavelength which can deviate slightly from £. This ansatz 
with Eqs. (2.3) and (3.5) imply [24] 

w^{x,t) = wo{l{x,t)), (3-8) 

i.e. only depends on the local fluctuations of cell size l{x,t), or alternatively on the 
cellular wavenumber, k[x^t) = 27r//(a;,t). Hence the time dependence of is given by 

where a is the response of the single-cell drift rate to a small change in cell size as given in 
Eq. (2.4). 
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The dynamics of k[x^t) can be obtained from symmetry considerations. For slight dis- 
tortions around the cellular solution, the wavenumber must be locally conserved. (Noncon- 
servation must involve the creation or annihilation of cells which is excluded in the cellular 
state by definition.) In addition, the inherent Galilean invariance of the KS equation must 
be respected. From these considerations, it follows that to leading order the equation of 
motion of k[x^t) must have the form 

dtk{x, t) = Xd^h^d^k{x, t) + \k{x, t)dlh< + /3a>< + -fdlk. (3.10) 

Note that the local conservation law for k excludes a term of the form —ak[x^ t) in Fq. (3.10). 
Comparing the second and third terms on the right hand side of Fq. (3.10), we find the 
constant fi to have the dimensions of A • A;. However, there is no preferred value of k[x^t) 
in the perturbation theory [25], so we must have /3 = 0. The constant 7 involves the cell 
height and is generically nonzero; its computation is more complicated and is not discussed 
here. 

To make the slow mode equations more transparent, consider A;(a;, t) = [2Tr /£)[l-\-6[x^t)]^ 
where 6[x^t) describes a small compression or dilation of the cellular solution of length £. 
Clearly, this new variable 6[x^t) is proportional to the local cell density. In terms of 8{x,t), 
the slow mode equations (3.4) and (3.10) becomes 

d,h^ = -vdlh^ + wo{l) + a{l)8{x, t) + [dth^,{d,h^f] , (3.11) 
d,8 = \dlh^+^dl8,+0 [djd,h^,8dlh^] , (3.12) 

where is the drift rate response. Although the above equations are derived in a phe- 
nomenological way, they can be obtained more systematically from a multiple scale analysis, 
as was done in Ref. [16]. 

Fquations (3.11) and (3.12) give a physical description of the slow mode dynamics. 
According to Fq. (3.12), cells slide into valleys {d^h^ > 0) and thereby become compressed 
{dt8 > 0), and similarly slide off of hilltops (9^/i< < 0) and become dilated {dt8 < 0). This 
behavior is a straightforward consequence of Galilean invariance as illustrated in Fig. 4. 
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Fluctuations in the local cell density affect the height profile through a fiuctuating local 
drift rate aS. For cr > in Fq. (3.11), cell compression in the valleys leads to a faster local 
drift rate pushing up the valleys, while cell dilation on the hilltops leads to a slower local 
drift, pushing down the hilltops. 

The above dynamics give rise to oscillations. If the terms with p and 7 (and all higher 
order ones) are neglected, then Fqs. (3.11) and (3.12) are easily combined to yield the wave 
equation 

d^h^ = c'dlh^, (3.13) 

where = Xcr[£). (It can be shown that the expression for given in Ref. [16] can be 
transformed to this form.) Thus for cr > 0, the cellular solution is elastic and supports wave 
propagation. From Fig. 3, we see that this occurs for £ ^ /*. For smaller £'s, < and the 
cellular solution is unstable to small perturbation. Thus for cr < 0, slower drifting smaller 
cells eventually become annihilated by their larger neighbors. 

In the regime > 0, the oscillation is weakly damped if 7 > leading to the behavior 
of visco-elasticity first found in Ref. [16]. However the oscillation becomes unstable if 7 < 
This was found to occur for £ > 0.91/o [16]. Thus stable visco-elasticity occurs for a narrow 
window of cell sizes, I* < £ < 0.92/o. 

We re-emphasize that stability of the large scale dynamics requires elasticity and damp- 
ing. Flasticity is governed by cr, which characterizes the response of the local drift rate to a 
small local compression. The analysis of this section establishes the (somewhat surprising) 
link between this response and the cr describing the size dependence of the single-cell drift 
rate. This is the case because changing the size of a single, periodic cell is an effective way of 
mimicking local compression in the cellular state of a large system. It is an example of how 
large scale dynamics can be predicted by computing some properties of a small, periodic 
system. We utilize this idea in Sec. IV for the chaotic regime. 

While an elastic response exists for a large range of parameters, damping is found to 
be weak even in the stable regime 7 > In particular, there is no diffusive relaxation as 
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expected of KPZ-type dynamics. This is a consequence of cell-conservation which excludes 
a DC-damping term of the form —ak in Eq. (3.10). Note that stability of the cellular state 
is possible if a damping term —fih is included in Eq. (1.1). Although such a term breaks 
the translational symmetry h ^ h -\- Const and is not allowed in the KS equation, it does 
describe a number of other interesting physical phenomena [26] and will be investigated 
elsewhere. Another possible route to stability is to break the A;-conservation. This requires 
cell creation and annihilation occurring in the chaotic state. 

IV. SPATIOTEMPORAL CHAOS 

The concepts developed for the perturbed dynamics of the cellular solutions provide a 
launching point for analyzing STC. As argued in Ref. [7], the system rarely ever settles 
into a cellular state if it starts from a generic set of initial conditions. Instead, it exhibits 
chaotic dynamics with inhomogeneous spatial structures as shown in Eig. 1. Chaos occurs 
for systems of sizes as small as 3/o ~ 4/o (see Refs. [14] and [15] for an in-depth discussion of 
the chaotic dynamics in small systems). The origin of chaotic motion in small systems has 
already been alluded to in Sec. II: The system is "frustrated" by the existence of multiple 
"metastable" cellular solutions. 

The chaos in a small KS system can be described in terms of its statistical properties 
given by the probability distribution ^[-ff], where H[x^t;£) is obtained from Eq. (1.1) in a 
periodic box of size £ with random initial conditions. In this section, we shall assume that 
VlH] is already known for £ of the order of a few /q's [14,15], and use it to obtain the statistical 
properties of the solutions in systems of sizes L ^ £. In Sec. V, we will numerically measure 
these statistical properties. It may be possible to attain some analytical understanding of 
the small system via amplitude equations [21]. The basic premise is that the STC of a 
large system can be understood in terms of the chaos observed in (mutually coupled) small 
systems. The size of the small system £ should be large enough to support a continuum of 
chaotic states. Let the shortest £ satisfying this criterion be Ic- We find in Sec. V that the 
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procedure worked well for £ > ~ 4/o. 

We wish to understand which properties of the small system are important, and how the 
small systems are coupled together. We accomplish this with a coarse-graining procedure 
similar to that used in Sec. Ill, dividing the /i-field into "slow" and "fast" parts, i.e. h = 
^< + with respect to £. The equation of motion for and hy are again Eqs. (3.4) and 
(3.6). It is important to note that the dynamics of the slow modes affect those of the fast 
modes through the last term in Eq. (3.6). 

To complete the description of slow mode dynamics, we need to specify the probability 
distribution of t/;<, obtained from the distribution of hy through the definition Eq. (3.5). 
As in the cellular case, we would like to relate hy to i7, the solution in a small, periodic 
system. However, generalization of the simple ansatz (3.7), 

hy{x,t) = H{x,t;l{x,t)), (4.1) 

no longer applies here, because a) "wavelength" is no longer well defined in the chaotic case, 
and b) the chaotic dynamics depend only weakly on the system size (provided the size is 
large enough to support the chaotic state). In particular, we can no longer mimic a local 
compression in a large system by changing the size of a small periodic system as was done 
in Sec. III. 

In the chaotic case, the simplest ansatz possible is 

Prob[/i>(x, t)] = V[n{x, t; £); d^h^] (4.2) 

where 7i is the solution of Eq. (3.6) in a small, periodic system of size £ under the infiuence 
of a slow imposed curvature d^h^. The explicit dependence of V on the slow mode comes 
from the last term of Eq. (3.6). (Since a constant d^h^ merely amounts to a Galilean boost, 
the leading nontrivial behavior comes from the dependence on d^h^.) Note that in using 
periodic boundary conditions for ?i, this ansatz neglects spatial couplings of hy between the 
neighboring coarse-grained regions. In this and the following sections, we will first derive the 
consequence of the ansatz (4.2), and then show self-consistently that the spatial couplings 
neglected do not alter the hydrodynamic behavior. 
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Given the ansatz (4.2), the probabihty distribution of becomes 

P[w<;dlh^] = Jvn6(w^- ^W<y^ nn; d^h^], (4.3) 

i.e. it becomes an exphcit functional of the chaotic solution of a small KS system, and can 
therefore be obtained numerically (see Sec. V) or possibly even analytically. Moments of 
the distribution P can then be used to construct an appropriate stochastic equation for w^. 
A review of this construction is given in the Appendix. Here, we shall motivate the form of 
the effective equation of motion for . 

From the definition of in Eq. (3.5), we have 

dtw^x.t) = x[d,hyd,{dthy)U, (4.4) 

where [•]< denotes spatial average as in Eq. (3.2). Substituting Eq. (3.6) for the dynamics 
of the fast modes leads to 

dtw^ = F[w^] + x^[d,hyd,{d,h^d,hy)U, (4.5) 

where all the purely fast mode contributions are represented by F[t/;<]. To find the form 
of F[t/;<], consider first the simple case where the system size is small, Z = £ ^ /c, so 
that the second term in Eq. (4.5) vanishes (i.e. no slow mode contributions). In this case, 
w^[t) = w^[x^t) is just the instantaneous, spatially-averaged drift rate of the small chaotic 
system, fiuctuating about its time-averaged value, {w^[t)) = vi. It is then natural to expect 
the dynamics of to be described by the following Langevin equation, 

dtw^{t) = -ai{w4t) - Vi) + C(t), (4.6) 

to leading order in — vi, with a response time and a Langevin noise £,{t). We will 
show in the Appendix that the numerical result of Sec. V leads to Eq. (4.6) as a correct 
description of the leading moments of P[t/;<]. The parameter a^, vi and the second moment 
of the noise 

mm) = Mt - 1') (4.7) 
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are determined numerically in Sec. V. 

The second term in Eq. (4.5) vanishes as well in a large system with L ^ Ic ii we set 
/i< = Const. In this case, we may expect the dynamics of to be similar to Eq. (4.6), but 
with a space-dependent and i^, i.e. 

dtw^{x,t) = F[w^] = -ai{w^{x,t) - vi) + C(a;,0- (4-8) 

Note that spatial coupling terms such as d^w^ should generically be present in Eq. (4.8). 
These are precisely the couplings neglected by the use of ansatz (4.2). They will be shown 
to be irrelevant in the hydrodynamic limit and will not be pursued here. In the simplest 
scenario, we expect the Langevin noise ^{x^t) to be uncorrelated in space for spatial scales 
> /c, i-e. 

{Cix, t)Cix', t')) = S{x - x')AfXt - f). (4.9) 

If this is true then we would expect the spatial average of i{x,t) (i.e. ^^{t)), to be given by 
the following correlator, 

{UrTiLm = Al(t) = Mr)/L, (4.10) 

for L exceeding several /c's. This l/Z-dependence will be verified by the numerics in Sec. V. 

Next we consider the effect of the second term on the right hand side of Eq. (4.5). It can 
be rewritten as 

x'iid^hyYdlh^U + Y^id.h^)dAd.hyyu 

= 2Xw^dlh^ + Xdxh^dxW^ + X[xdxw]^dlh^. (4-11) 

The dominant term is obtained by approximating by its time average vi. Then to lowest 
order Eq. (4.11) becomes 2Xvfd^h^. Higher order terms are again irrelevant as will be shown. 
Putting the above together, and re-expressing in terms of the local drift rate fiuctuation, 

TT^X^t) = w^{x,t) — Vl, (4-12) 
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we obtain the equations 

d,h^= -vdlh^ + ^{d^h^f + vt + 7r(x, t) + O [a><] (4.13) 

= -atT + fitdlh^ + i{x,t) + o\dlT,d,Td,h^,Tdlh^,T\{dlh^f] , (4.14) 

describing the coarse-grained dynamics, with fii = 2Xvi. A formal construction of Eq. (4.14) 
based on the numerical results of Sec. V is presented in the Appendix. 

It is useful to compare Eqs. (4.13) and (4.14) to the cellular slow mode equations (3.11) 
and (3.12). We note that the drift rate fluctuation Tr[x^t) resembles the cell density fluc- 
tuation 8{x,t). This should be expected: For any quasi-cellular profile the local drift 
rate operator w^{x^t) = ^^{dxhyY^ has two peaks in x for each cell, therefore the long 
wavelength part of w^[x^t) should provide a measure of the local cell density (at least in 
one spatial dimension). Identifying Tr[x^t) with the cell density fiuctuation, we see that as in 
the cellular case, the chaotic dynamics contains elasticity. Cells are compressed [dfTr > 0) in 
valleys {d^h^ > 0) and dilated [dfTr < 0) on hilltops {d^h^ < 0). Compressed (dilated) re- 
gions have higher (lower) drift rates, thus providing a restoring force to height fiuctuations. 
The elasticity is again quantified by the drift response of the system to compression. In 
the cellular case, compression was effectively mimicked by changing the cell size, but in the 
chaotic case, compression is produced by explicitly imposing an external curvature d^h^ as 
described in Sec. V. (Defining aS = tt^ we see that the drift response a = 2vi is always pos- 
itive in the chaotic case.) The main additions due to the chaotic dynamics are the damping 
term — a^vr and the noise term ^{x^t) in Eq. (4.14). As alluded to earlier, the damping term 
is allowed in the chaotic dynamics because cell density is no longer a conserved quantity 
(there is a preferred A;), while the noise ^{x^t) describes the stochastic nature of the cell 
creation and annihilation process. In fact, Eq. (4.14) is the minimal equation demanded by 
symmetry in the chaotic regime. 

The damping term makes the chaotic dynamics fundamentally different from the cellular 
dynamics. In the hydrodynamic limit, Eq. (4.14) can be solved to yield 
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(4.15) 



Substituting Eq. (4.15) into the slow mode equation (4.13), we finally obtain the effective 
equation 

dth^ = v.dlh^ + \{d,h^f + v, + r,{x,t) + o [dth^,dld,h^, [dlh^Y] , (4.i6) 

with 

l^i = ^-i^. (4.17) 

{ri{x,t)ri{x',t')) = S{x - x')Ai{t - t')/al (4.18) 

As will be shown in Sec. V, A^(t) is correlated only up to a time scale Tq ~ to. Thus for time 
scales ^ To (in the hydrodynamic limit), it may be approximated by a (5-function. Hence 

{riix, t)T]{x\ t')) = Df8{x - x')8{t - f), (4.19) 

with an effective noise amplitude 

= / dTAi{T)/al (4.20) 

The "bare" parameters Vf ^Df , are determined in Sec. V through simulations of small KS 
systems of size L We find the result Di > which establishes the hydrodynamic system 
Eqs. (4.16) and (4.19) to be in the KPZ universality class, provided that the coefficients of 
the irrelevant terms (those in the square bracket of Eq. (4.16)) are not pathological. In order 
to compare the results to the parameters z/gff, -Defr of the effective KPZ equation measured 
at very large scales [12], additional renormalization of vi^Df by the irrelevant terms must 
be considered. This issue will be addressed in Sec. VI. 



V. NUMERICAL METHODS AND RESULTS 

The effective hydrodynamic equation given by Eq. (4.16) and Eq. (4.19) must be checked 
for self consistency by computing the parameters vi and Df . One must demonstrate that the 
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effective surface tension vi is positive and tliat tlie effective noise is spatiaiiy and temporaiiy 
uncorreiated. This is accompiished by malcing numericai measurements of tlie ieading mo- 
ments of tlie probability distribution function P\w^] 9j/i<]. As shown in the Appendix and 
in Sec. IV the moments are directly related to the parameters of the effective equation. The 
task is greatly simplified by the use of ansatz (4.2) which replaces a small segment of a large 
system by a small periodic system of the same size subject to some relevant external forcing 
(curvature). The spatial couplings neglected are recognized as irrelevant once Eq. (4.16) is 
established to described the hydrodynamic limit. 

In order to compare consistently to the results of the large KS system obtained in 
Ref. [12], we used the same numerical scheme: A simple real space finite differencing 
scheme was used with the nonlinearity discretized as |-[(/ii_|_i(t) — /ii_i(t))/2]^, where i is 
the discretized spatial index. An Euler step was used for the time evolution, with step size 
di = 0.1. The simulations were done with the same parameters as those in Ref. [12], namely 
;/ = 1, /€ = 1, A = 2, for which /q ~ 8.89 and to = 4. The measurements were made on 
systems ranging in sizes from L = 32 to L = 150, corresponding to 3.6/o to 16.9/o. Simula- 
tions were performed on an IBM RS6000 workstation unless otherwise specified. It should 
be noted that the lattice KS equation studied in Ref. [12] and here is actually a dynamical 
system distinct from the continuum KS equation (1.1). Although we find that the lattice 
KS equation belongs to the same universality class as the continuum equation, eventually 
fiowing to the same KPZ fixed point, the precise values of the "bare" parameters of the 
resulting effective hydrodynamic equation (4.16) do depend on numerical and discretization 
schemes. However, since our approach is not very computationally demanding, a careful, 
high accuracy numerical method could be easily employed to get values that better represent 
the continuum KS equation. 

As shown in the Appendix, it is necessary to obtain the first and second moments of 
the distribution P[t/;<; 9j/i<], as well as the response function. In addition, due to the 
discretization scheme used, Galilean invariance is broken at the microscopic scale. This 
leads to a modified response to uniform tilt d^h^ (a renormalized nonlinear parameter Agfr) 
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which is also investigated. 



1. Response Function and Stability 

The effective surface tension may be obtained Irom the time-dependent response ol the 
dnlt rate to an imposed curvature d^h^ in small systems of the size of several /q's (see 
Appendix). Recall from Galilean invariance that a slow curvature compresses or dilates the 
cells (see Fig. 4), thereby changing the drift rate. It is this change of drift rate in the chaotic 
state that is desired. The measurement was implemented numerically by simulating the 
dynamics of the fast variable Eq. (3.6) with a slow curvature imposed through the choice of 

Ideally a uniform forcing d^h^ = Const (i.e. oc x"^) should be applied. However, in 
practice an infinitely large parabola can never be used and boundary effects will become 
important. The drift response will not be uniform at the boundaries thus distorting mea- 
surements. This problem may be overcome in a number of ways. One could use a large 
enough parabola and only look at the drift response near the center away from the bound- 
aries. However, depending on the length scale over which the response is probed, this could 
require a rather large system and become computationally intensive. Another way is to ap- 
ply a slowly varying d^h^ (e-g- ^ sine wave), systematically remove the nonuniform effects 
and measure the response to the average curvature imposed. We choose the second method. 

With a spatially nonuniform forcing Eqs. (4.2) and (4.3) must be slightly modified. 
One should use = ^[d^'Hy Y^ where Tiy comprise of the Fourier modes that are "fast" 
compared to the spatial variation in d^h^. The corresponding "slow" modes ?i< are a 
"back-reaction" to the nonuniform curvature and must be subtracted away. The back- 
reaction arises because each region experiencing a different curvature will have a different 
response. This will then distort the surface and alter the effect of the imposed curvature. To 
ensure that the back reaction has been properly removed the time averaged spatial profile 
must be monitored. The criterion for effective back reaction neutralization is that the time 
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average of the spatial profile remains flat. i.e. {dxTi-{x^t)) = 0. 

A sine wave was specifically chosen for the slow forcing because the shape of the "back 
reaction" will also be a sine wave to linear order and can be easily removed. However, there 
are still nonlinear effects that can alter the background curvature. If significant spatial 
nonuniformity in {7i) still exists after the first Fourier mode is removed, then higher modes 
are also removed until the surface is uniform to an acceptable tolerance (< 10% of the forcing 
amplitude). For the system sizes used we never needed to remove more than two modes. 

In the simulation, Fq. (3.6) was used with 



Hence the left and right halves of the system were under the infiuence of a positive and a 
negative curvature respectively. In each run, two measurements of the drift rate, one for 
the left and one for the right half were made. We took the imposed curvature c to be the 
'average' curvature experienced by each half. 



The first Fourier mode of /i>, hy[k = vr/Z), was removed after each time step. This was 
accomplished by using an FFT, zeroing the first mode, then using an inverse FFT to restore 
to real space. The DC mode was also zeroed to remove the average drift component. The 
second mode was removed if necessary. 

The measurement was performed by starting the simulation with c' = and with small 
random initial conditions. After the system settled into the chaotic attractor (several linear 
growth times), a small curvature (c' > 0) was imposed at t = and the new drift rate was 
measured by taking w^[t) = ^[j^J^idxhyY] separately for the left and right halves. The 
measurement was then averaged over many configurations of random initial conditions for 
each used. Our main measurement was for Z = 64 averaged over 10^ samples and for c' 
ranging from to 0.004. This corresponded to Z ~ 7/o which was large enough to capture 
the essential small system chaotic behavior yet small enough to be computed comfortably. 




for 
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The spatial nonuniformity after the removal of the first Fourier modes was found to be under 
10% of the amplitude of h^. We also made measurements for L = 32 and L = 128 but with 
fewer samples and a smaller range of c'. 

Fig. 5 shows a typical form of {w^[t)) for c' = 0.002 and Z = 64 for the two halves of the 
system. The response rises/drops sharply before eventually saturating to a constant value. 
This behavior is well described by the exponential form 



to be independent of c', while A(c') increases monotonically as shown in Fig. 6. A least- 
square fit finds 



The effect of the nonlinear term is small over the range of c' used. From the exponential 
form of the linear response Fq. (5.3) we obtain (see Appendix) the effective equation of 
motion for , 




(5.3) 



The values of Al[c') and aL(c') were computed. We find 



aL = 0.15 ±0.02 



(5.4) 




(5.5) 




(5.6) 



where c = d^h^ is the average slow curvature forcing. We can identify 




(5.7) 



which leads to 





using c' = TTC from Fq. (5.2). This gives an effective surface tension of 



Vl = — -v = 13.5 ±0.6. 



(5.9) 
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These measurements have been repeated for different Z's. For L = 32 we found the 
values: a/, = 0.15ib0.03, z/L = 12.5ibl. For L = 128, removing the first Fourier mode was 
not enough to counter the back reaction. We removed the first two modes and obtained the 
results: a = 0.14 it 0.03 and Vl = lA i: 2. The larger uncertainties are due to the fact that 
fewer samples and a smaller range of c' was used in the measurements. We find only a weak 
Z-dependence for the parameters and jS^ as long as the small KS system is well within 
the chaotic state and the back reaction is removed. 

We have also tried a different forcing term consisted of two parabolas of opposite curva- 
tures: 



instead of the sine forcing Fq. (5.1). In this case the back reaction was removed by removing 
the first mode and/or removing the average drift of the right and left halves. We find results 
similar to those reported above, except with larger uncertainties. This is attributed to the 
fact that our procedure of removing the back reaction was incomplete even to linear order. 



The effective noise amplitude was calculated from the two-point correlation function of 
the drift rate t/;<(t). 



for a small system with no external forcing and with periodic boundary conditions. The 
numerical scheme used was the same as that for the response function. Starting from random 
initial conditions and after allowing the system to reach the steady state, Cl{t) was measured 
up to a time when the correlation function saturates. The result was then averaged over 
10"* configurations of random initial conditions. A typical form of the averaged correlation 
function is shown in Fig. 7 for L = 100. From Fq. (5.6) with c = and for intrinsic noise 




for < x < I 
for ^ < X < L 



(5.10) 



2. Noise Amplitude 
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(^(x, t) that is completely uncorrelated in space and time, ((^(x, t)^[x' ^ t')) = B8{x — x')8{t—t'), 
then 

The fine structure in the numerically obtained Cl{t) shown in Fig. 7 indicates short-time 
structure in (^(t), which can be obtained by de-convolving Cl{t) from the linear response 
function G{t) as in Eq. (A18). However, this short-time structure does not affect the 
hydrodynamic limit, and is not pursued. As explained in Sec. IV, the quantity relevant to 
hydrodynamic behavior is the time integral of the correlator 

Dl{t) = {{w4t) - VL){w<{t + r) - vl)) = i (C4^) - Cl{t)) . (5.13) 
If is short-range correlated in space, we expect 

/'CO 

Dl = L dTDLir) (5.14) 

J — CO 

to be independent of L. This quantity has been measured for various system sizes ranging 
from 50 to 150. We find 

1)^ = 2.6 ±0.1, (5.15) 
with negligible Z-dependence (see Fig. 8). 

3. Breaking of Galilean Symmetry 

In Ref. [12] it was found that the breaking of Galilean symmetry by finite differencing 
resulted in an increase of A from the bare value of 2 to 4.69 it 0.15 for a large KS system 
with L = 1024. We performed the same measurement by determining the dependence of 
the average drift rate vl to a macroscopic tilt 6, which is imposed on the system via helical 
boundary conditions, h{L) = h[0) + 6L. Fitting the data to the parabola vl{0) = + 
we found \l = 4.4 ± 0.2 and vl = 0.493 ± 0.001 for L = 100. The drift velocity vl has a 
very weak dependence on L. As in Ref. [12], we found A to be sensitive to both the time and 
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spatial discretization size, with Xl approaching A as the discretization step sizes are made 
smaller. 

Another way of determining the effect of Galilean symmetry breaking is through a mea- 
surement of the advection rate in the presence of a macroscopic tilt. For the continuum 
system, Galilean invariance gives x ^ x -\- ut, with the advection rate u = \6 ii h ^ h -\- 6x. 
However, since our discrete system breaks Galilean invariance, it is not a pn'on obvious what 
u{6) should be in this case. In the chaotic state, one can measure u{6) through the average 
response function {Q{x,t)), which can be obtained by following the difference in time evolu- 
tion of two samples with identical (generic) initial conditions except for a small difference at 
the origin (x = 0), and then averaging over all pairs of random initial conditions. For ^ = 0, 
{Q{x,t)) is an oscillating function in x, peaked at x = 0, with decaying amplitude for large 
X (see Fig. 9). With ^ 0, the advection rate can be measured through the movement of 
the peak of {Q[x^t)) in x. In practice, however, the above is a very time-consuming task, as 
a large ensemble average is needed to obtain {Q[x^t)) for each 6. 

To get an idea of the behavior of the advection rate, we measured the tilt-dependence 
of our discretized system for the stationary cellular solution. First we checked the tilt- 
dependence of the drift rate as defined in the previous paragraph. We found A = 2.8 it 0.2 
for L = 8 and A = 4.8 it 0.2 for L = 10. This shows a strong lattice effect in the cellular 
solutions. (The values of A obtained in the stationary and chaotic states are not expected 
to coincide, since as already described in Sec. Ill, the drift rate of the stationary state has 
a strong Z-dependence even for 6 = 0.) Nevertheless, when we measured the advection rate 
we recovered the result u = X6 with the bare value of A independent of L. We believe this 
result concerning the non-renormalization of the advection rate to hold also in the chaotic 
state, since the local advection rate there should be similar to the size-independent advection 
rate of the stationary state. 
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VI. ANALYSIS OF RESULTS 



The existence of a linear response of the drift rate to external curvature d^h^ (Fig. 6) 
validates Eq. (5.6) as the equation of motion for t/;<, or equivalently Eq. (4.14) for the drift 
rate fluctuation Tr[x^t) = w^[x^t) — v^. This leads to Eq. (4.16) as the effective equation of 
motion for the slow modes h^[x^t)^ with the effective surface tension 



which depends weakly on £. The positivity of Di along with the independence of Di on £ 
(Eig. 8) establishes Yakhot's conjecture that the KS dynamics belongs to the universality 
class of the KPZ equation. It then follows that the terms neglected in our treatment, in 
particular, the spatial coupling terms (Sjvr, etc. in Eq. (4.14)) neglected due to the use 
of the ansatz (4.2), are irrelevant at the KPZ fixed point. We of course cannot a priori 
exclude the possibility that some of the coefficients of the irrelevant terms are very large, 
e.g. the neighboring chaotic attractors are strongly coupled spatially, thus driving the system 
to a different fixed point. This possibility can be checked by applying a space-dependent 
curvature perturbation in Eq. (3.6) and measuring the space-dependent response of the drift 
rate. The lack of strong ^-dependence in the parameters Vi and Di measured suggests that 
such effects are not important. 

The remainder of this section is devoted to a more quantitative analysis of the numerical 
results obtained. Erom the measurement of the saturated drift rate response to external 
curvature (Eq. (5.5)), and the response time ai (Eq. (5.4)), we obtained jSi = 2.2 it 0.2. The 
analysis of Sec. IV predicts 



The drift rate for for £ = 64 is 0.49. Thus Eq. (6.2) is satisfied within numerical 

uncertainties if the bare value A = 2 is used. There is, of course, an ambiguity as to 
which is the "proper" A to use for the lattice KS equation studied here, since the effective 



Pi = p ^ 



13.5 ±0.6. 




/3i = 2Xvi. 
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A obtained from drift and advection response to tilt are different as shown in Sec. V. We 
argue that the proper A to use is the advection response, A = 2. As already explained, 
the drift rate fluctuation tt is proportional to cell density fluctuation which responds to a 
nonzero background curvature through advection. The cell density increases at the bottom 
of a valley because the advection process brings cells in from the two sides of the valley (see 
Fig. 4). Thus, the response of cell density to changes in background curvature is determined 
by the advection rate, which is proportional to the bare A. (A similar ambiguity also exists 
for perturbations about the cellular solutions, where the analysis based on the continuum 
equation predicts wave dynamics Eq. (3.f3) with (? = Act. For the lattice KS equation 
studied here [and possibly the pseudo-spectral system studied in Ref. [f6] although A was 
never measured], it is again unclear as to which is the proper A to be used. However, our 
numerics [and those of Ref. [f 6]] clearly support the use of the bare A and are incompatible 
with the renormalized value A obtained from the drift response.) 

As mentioned in Sec. Ill, the large scale KS equation is parameter free in terms of the 
"natural" units, the typical cell size /q = 27r^2/€/z/ 8.89, the linear growth time (which 
is of the order of the maximum Lyapunov time) to = Ak./p'^ = 4, and the cell height ho = 
2v I \ = I. In these units, the bare parameters are v 0.7 [/q/^o], D 1.2 [K^lojio], i.e. they 
are of order unity. This should be expected from the defect-mediated stability mechanism — 
both the stability factor v and noise D are due to the microscopic cell creation/annihilation 
process which involves the change of the local interface height by an amount ho over a spatial 
extent of /q and is completed within a time to- 

For the same reason, we expect the bare coefficients of the irrelevant terms, (terms in 
the square brackets in Fq. (4.16)), also to be of the order unity in their respective natural 
units, at the smallest possible scale Ic at which Fq. (4.16) can be defined. These irrelevant 
terms will certainly cause finite renormalization of the effective parameters vi and Df (by 
an amount of the order unity in natural units), eventually leading to the effective KPZ 
parameters z/gff and Dgff measured in Ref. [12] for coarse-graining scales £ ^ /c. Thus the 
parameters vi = 13.5 it 0.6 and Df = 2.6 it O.I measured in Sec. V are not quite the 
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macroscopic parameters z/gff and Dgff since £ = 64 used is only about twice Ic- They serve as 
an order-of-magnitude estimate of z/gff and -Deff- Better estimates of the effective parameters 
can possibly be made by repeating the numerics of Sec. V for larger systems. However, since 
renormalization by the irrelevant terms quickly saturates for £ > Ic, the weak ^-dependence 
of the parameters vi and Df found suggests that the range of £ used is probably already 
close to the saturated region. Thus the numerical values of vi and Df obtained should be 
close to z/gff and -Deff- This is indeed the case, since Ref. [12] finds z/gff = 10.5 it 0.6 and 
i)eff = 3.2 ±0.1. 

It should be noted that an accurate determination of the bare parameters of the KPZ 
equation is possible in this case only because the cross-over length scale beyond which 
the KPZ nonlinearity becomes important is unusually large. As estimated in Ref. [12], 

/x ^ 150/o/?, (6.3) 



where the dimensionless coupling constant is 



'^eff 



/o ^ 0.53. (6.4) 



Since /x ^ ^0 ~ ^c, it is possible to perform the numerics at £ ^ /c, and include the 
effect of the irrelevant operators without running into the complication of renormalization 
by the relevant KPZ nonlinearity. Note that the large value of is a result of the large 
prefactor, 150 in Eq. (6.3). This is purely a "KPZ effect". What refiects the inherent nature 
of the KS dynamics is the dimensionless constant g which is of order unity. Of course, 
a large separation of scales between and /q should not be expected to hold for generic 
STC problems. Without a clear separation of scales, it will not be possible to determine 
the "bare" parameters of the effective hydrodynamic theory by this or any other methods. 
(Indeed, the bare parameters will not even be observable in such cases.) Nevertheless, one 
can still use the approach given here to construct the form of the effective equation, from 
which the asymptotic, universal behavior may be obtained. 
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VII. SUMMARY AND FUTURE OUTLOOK 



We have presented an effective hydrodynamic theory of spatiotemporai chaos exhibited 
by the Kuramoto-Sivashinslcy equation. We constructed an effective stochastic equation 
beionging to the KPZ universaiity ciass in the hydrodynamic iimit by incorporating the 
chaotic dynamics of the smaii KS system in a coarse-graining procedure. This resuit was 
estabiished by computing the bare parameters of the effective equation. They are found 
to be of the same order as the ones obtained from large scale simulations. An important 
assumption in our approach is weak spatial coupling of the local chaotic dynamics. This 
allows the use periodic boundary conditions when measuring properties of the small KS 
systems. The assumption is justified because of the irrelevancy of these spatial couplings at 
the KPZ fixed point. 

Due to Galilean invariance, fiuctuations in the background curvature are coupled to 
fiuctuations in the cell density. Thus a preference for zero curvature (positive surface tension) 
is unavoidably connected to a preference for a specific cell density. The latter is only possible 
in the presence of spacetime defects, which are manifestations of the division of thick cells 
and the coalescence of thin cells. We refer to this process as "defect-mediated stability". 
The dynamics of the defects, which involve the local change of wave number, are chaotic 
in the KS system because a band of states with nearby wave numbers are metastable and 
thereby frustrated as discussed in Sec. II. Thus stability is accompanied by stochasticity. In 
the absence of defects (the cellular case) wavenumber conservation only leads to visco-elastic 
motion at best. 

It is instructive to examine KS-like equations, i.e. equations similar to Eq. (1.1) but with 
linear operators of the form — b > a. From symmetry, one might again conjecture 
that these systems belong to the KPZ universality class. However, in light of the analysis 
presented, the stability mechanism can only produce a term proportional to —k^ in the 
slow mode equation. Therefore, one can conclude right away that a necessary condition for 
stability is a > 2. This explains why instabilities from the linear operator \k\ — cannot 
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be tamed by (dxh) -type nonlinearity, a result well known from the study of the viscous- 
fingering instability [27]. On the other hand, KPZ behavior is expected for a > 2. The bare 
coefficients can be computed using the scheme of Sec. V and can be compared to those from 
large scale simulations [28]. 

If the system is such that the band of stable, single-cell solutions (the analog of the range 
of sizes between li and I2 in Sec. II) is very narrow compared to the typical cell sizes, then 
the size of the "small" system required to extract local chaotic dynamics must be very large. 
This is because frustrations due to the coexistence of multi-cellular solutions of different 
wavenumbers is needed for chaos (recall the four-cell, five-cell example in Sec. II), while 
the minimum system size Ic required to accomodate different wavenumber states is inversely 
proportional to the stable band width. It is tempting to conjecture that the density of 
positive Lyapunov exponents will scale as 1/lc at length scales ^ Ic- Such behavior may be 
investigated in systems exhibiting an order-chaos transition. 

It should also be interesting to extend the present approach to study the behavior of 
the 2-|-I dimensional generalization of the KS equation. Again from symmetry, one can 
make a case for the KPZ equation to describe the hydrodynamic behavior. This approach is 
advocated by Jayaprakash et al. [29]. On the other hand, fundamentally different behavior 
is proposed by L'vov and Procaccia [30]. The differences between the two proposals are 
not expected to be resolved numerically with current computational capabilities. However, 
it may be possible to construct the effective hydrodynamic equation explicitly with our 
method. A straightforward generalization of the ID analysis points to some important 
differences between the ID and 2D systems: The configurations of the cells can no longer 
be described by a scalar density field, since the local orientational degrees of freedom must 
be specified in 2D. Also, the local cell arrangements are determined by the local curvature 
tensor, leading to more complex elastic responses. The minimal effective equation one can 
construct in this case is the analog of Eqs. (4.13) and (4.14), with Tr[x^t) generalized to 
a complex field, whose phase describes the orientational degrees of freedom. Recall that 
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symmetry is valuable in constructing the effective hydrodynamic equation only if all of the 
slow degrees of freedom have been identified. The existence (and possible relevance) of the 
orientational degrees of freedom in the 2D KS system therefore makes the correspondence to 
KPZ dynamics less transparent. A study of the 2D KS dynamics will be reported elsewhere. 

On a broader scope, it is important to see to what extent the approach presented can 
be applied to characterizing and understanding different types of STC. It seems most suited 
to describing phase turbulence where space-time dislocations or defects are interspersed in 
a background pattern that is otherwise periodic. Our analysis suggests that effective hy- 
drodynamic equations may be obtained by simply breaking wavenumber conservation and 
including a noise source in the appropriate slow mode equations describing the long wave- 
length distortions of the periodic patterns (the analog of Equations (3.11) and (3.12)). On 
the other hand, due to the weak spatial-coupling assumption used, this approach probably 
may not be helpful in understanding STC in coupled chaotic maps which rely crucially on 
the spatial coupling among neighboring maps. For such problems, the method recently intro- 
duced by Hansel and Sompolinsky [31] may be more appropriate. Through the exploration of 
these issues, one can hopefully gain more insight into the classification and characterization 
of various types of spatiotemporal chaos. 
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APPENDIX A: CONSTRUCTION OF THE LANGEVIN EQUATION 

We review the formal construction of a stochastic equation for a time series w[t) subject 
to a time-dependent perturbation e(t) given the leading moments of the distribution P[w; e]. 
It is convenient to consider the Fourier transform of P, 

P [td; e]= [Vw P[w; e] e"''/'^* (Al) 



which generates moments of P through the Green functions 

(5e (ti) • • • (5e bw {t.ra+i) ■■■bw 

For example, 

G^r{t) = {w{t))=v, (A3) 
^02(^1,^2) = {{w{h) - v) {w{h) - v)) (A4) 

are just the mean and variance of the distribution P[t/;; 0], and 

Giiih.h) = j^AHh)) (A5) 

be[ti) 

is the linear response function. Writing P = e"^, we see that Gmn are just the coefficients of 
the Taylor expansion of S: 

S [w, e] = „_o / dti ■ ■ ■ dtm+n j rGmnitlj ' ' ' 7 tm+n) 

mini 

e(^i) • • • e{t^)w{t^+i) ■ ■ ■ w{t^+n)- (A6) 

The "dynamic action" S can now be used to construct the effective equation of motion. 
Suppose the leading moments of P[t/;; e] [Fqs. (A3), (A4) and (A5)] exist, and the higher 
order moments are negligible. In this case, the action becomes 

S[w,e] = J dtdt' Go2{t, t')w{t)w{t') - iv J dt w{t) - i J dtdt' G'ii(t, t')e{t)w{t'). (A7) 

Rewriting the first term of Fq. (A7) by another Fourier transform, we have 
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P[w;e] = J Vt] exp|-z J dt [v + T]{t)]w{t) - i J dtdt' Gii{t,t')e{t)w{t')^ X P [t]] , (A8) 



where 

PM oc exp|-i J dtdt' G-Ut,t')r,{t)r,{t')^ , (A9) 

with Go2^(t, t') denoting the inverse of Go2{t, t'). Inverse Fourier transforming Eq. (A8) leads 
to the result 

P [w- e]= j VtdP [td, e] e'l^' '^W^'W 

= j Vri\{8 w{t)-v -ri{t)- j dt'Gii{t\t)e{t') X P [77] . (AlO) 

Eq. (AlO) clearly describes a time series which satisfies the "equation of motion" 

j dt' G^^ (t, f) [w{f) -v] = i{t) + e(t), (All) 

where 

i{t)= j dt'G-l{t,t')r,{t') (A12) 



plays the role of a Langevin noise. By the distribution (A9), we see that is Gaussian 
distributed, with the variance 



{({tW)) = Mt -t') = J dTdT'G-l{t,T)G^l{t' - t')G'o2(t - r'). (A13) 

In the presence of higher moments of P[t/;;e], it is straightforward to verify that P[vi\ = 
P[v-\-ri] 0]. In addition, the equation of motion (All) also acquires nonlinear terms. However 
for the KS equation, the nonlinear terms as well as the higher moments of P[vi\ are irrelevant 
in the hydrodynamic limit and are not pursued here. 

In Sec. V, we find the linear response of the average drift rate {wit)) of a small system 
to a curvature perturbation eit) = d^h^Oit) to be well described by the exponential form 

{wit)) = V + ^ [l - e-"-') dlh^. (A14) 

Since {w[t)) = v -\- J di'G\\{t\t)eit') according Eq. (AlO), we find 
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G^^{t\ t) = /3e-«(*-*')^(t - t'). (A15) 
Inverting Eq. (A15) yields the operator 

G-i (t, t') = 6{t - f) [dt + a) (A16) 
Finally, applying Eq. (A16) in Eq. (All), we obtain the equation of motion 

d^w = -a{w{t)-v) + I3dlh<+i{t), (A17) 

with 

mm) ^ A(t - t') = {a' - df)Go2{t - f ) (A18) 
to leading order in w[t) — v and i^. 
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FIGURES 
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FIG. 1. Spatiotemporal chaos produced by the KS dynamics. 
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FIG. 2. Shape of the cellular solutions for L/Iq = 0.72, 0.85, 1.0 (solid line), and 1.2. Both the 
width and height of the cells are normalized by L. 
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FIG. 3. Size-dependence of the single-cell drift rate wq{L). Iq = 8.89 is the most unstable 
wavelength according to the linear stability analysis, li = 0.707/o is the length where the nontrivial 
uni-cellular solution first appears. I2 = 1.27/o is the maximum size of a stable uni-cellular solution. 
wq{L) has a maximum si L = I* = 0.85/o. 



37 



FIG. 4. Galilean invariance requires the cells to advect in the downhill direction, leading to 
accumulation of cells in valleys and depletion of cells on hilltops. 



38 



0.500 



0.495 - 



A 0.490 - 



V 



0.485 - 




0.480 - 



10 



15 
t 



20 



25 



30 



FIG. 5. Time evolution of the average drift rate {w^(t)) after a small curvature forcing 
(c' = 0.002) is turned on at i = after the transients have died away. The data is taken for 
a 128-site system and averaged over 10^ configurations of random initial conditions. The upper 
(lower) branch corresponds to the drift rate measured from the left (right) half of the system. 
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FIG. 6. The saturated response Al(c') of an i = 64 system for various imposed curvature c'. 
The solid line is a least-square fit (see the text). 
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FIG. 7. The auto-correlation function Cl{t) of the drift rate of an i = 100 system. The data 
shown is averaged over IC* configurations of random initial conditions. 
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FIG. 8. The noise amplitude Dl as a function of L. The solid line shows no L dependence. 
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FIG. 9. The response of the KS equation to an infinitesimal perturbation of the initial condi- 
tions at a; = 0. The result is averaged over 51,200 configurations of random initial conditions. The 
measurement was performed on a CM2 connection machine. 
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